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PACS. 61.30.Pq - Microconfined liquid crystals: droplets, cylinders, randomly confined liquid 

crystals, polymer dispersed liquid crystals, and porous systems . 
PACS. 61.30.Vx - Polymer liquid crystals . 
PACS. 61.30. Jf - Defects in liquid crystals . 

Abstract. - We investigate condensation of a long confined chiral nematic polymer inside 
a spherical enclosure, mimicking condensation of DNA inside a viral capsid. The Landau-de 
Gennes nematic free energy Ansatz appropriate for nematic polymers allows us to study the 
condensation process in detail with different boundary conditions at the enclosing wall that 
simulate repulsive and attractive polymer-surface interactions. Increasing the chirality, we 
observe a transformation of the toroidal condensate into a closed surface with an increasing 
genus, akin to the ordered domain formation observed in cryo-microscopy of bacteriophages. 



Introduction. - DNA undergoes a pronounced change in volume upon addition of various 
condensing agents such as polyvalent cations (Spd^+, Spm^+, CoHex, polylysine, and histone 
proteins) but also in the presence of monovalent cations on addition of crowding agents such as 
polyethylene glycol (PEG) 1 . Under restricted conditions of very low DNA concentrations in 
the bulk, condensation proceeds via a toroid formation 2 that has been studied extensively 
with cryo-electron microscopy methods in order to determine its detailed morphology and 
internal DNA ordering 3 . Recently, toroidal condensates of a single DNA chain were observed 
inside an intact viral capsid, on addition of condensing agents Spm^+ or PEG to the outside 
bathing solution of the capsids Surprisingly, the shape of the condensate depended on 
the condensing agent, so that Spm^+ condensed aggregates looked very similar to the bulk 
toroids, whereas PEG condensed aggregates were flattened and non-convex, adhering to the 
capsid inner surface [5] . It thus appears that DNA-capsid wall interactions play an important 
role in modifying the morphology of a capsid-enclosed DNA aggregate, a statement that we 
will analyse in more detail below. 

Previously, the shape of DNA toroid confined to a spherical capsid shell was analyzed using 
the Ubbink-Odijk theory [7||8], that need to be modified appropriately in the case of attrac- 
tive DNA-capsid surface interactions [s]. All these theories, based on the elastic deformation 
energy Ansatz, treat the inhomogeneous DNA ordering only approximately and can not be ap- 
plied at all to describe the nematic transition itself. Nevertheless, this approach is completely 
appropriate to derive the shape of the already ordered DNA phase. Recently we proposed 
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a change in perspective and analyzed the confined nematic polymer such as DNA in terms 
of nematic ordering framework by writing down a Landau-de Gennes type confined nematic 
polymer free energy |10 . Alternative formulations of confined nematic polymer ordering can 
be based on density functional theory as in Ref. 11 or on a minimal, coarse-grained elastic 
model of densely packed confined polymer chains as in Ref. 12 . Motivated primarily by the 
recent experiments on condensed DNA states inside bacteriophage capsids 5,9, we gener- 
alize the confined nematic polymer analysis in three important aspects. First by explicitly 
investigating the condensation transition, then by adding the condensate - confining surface 
short range interactions that can be either repulsive or attractive, and finally by allowing 
for a chiral coupling in the free energy. We then study the dependence of the equilibrium 
order of the confined polymer on its density as well as the effect of short range attractions 
between the polymer and the confining surface on the nature of ordering and morphology of 
the condensate. 

Theoretical model - To determine the equilibrium configuration of the director and 
density fields we set up a free energy density following the Landau-de Gennes approach p!Qp3] . 
The appropriate variables in the polymer case are the complete non-unit nematic director field 
a(r), describing the orientation and the degree of order, and the polymer density field p(r), 
expressed as the volume density of chain segments of length £q. Both fields are coupled by 
the continuity requirement for the "polymer current" j(r) [10^: 

V-j(r) = p±(r), j(r)=/,(r)4a(r), (1) 

where is the volume number density of beginnings (p^ > 0) and ends (p^ < 0) of the chains. 
The conservation of polymer mass will be satisfied globally by requiring J dVp = mo = const. 
The ordering transition will be controlled by the density (concentration) of the polymer as 
is usually the case for lyotropic nematic liquid crystals [13|. For computational reasons all 
equations must remain regular for vanishing order, i.e., they must be expressed by the full 
vector a. Decomposition of the form a = an, where a is the degree of order, would result in a 
singularity of the form x oo in the centers of defects, where Vn diverges while the degree of 
order vanishes. In contrast, a and its derivatives remain regular everywhere. In this respect, 
Eq. ([T]) is already of the correct form. 

In the elastic free energy, instead of using the usual Frank terms for splay, twist, and bend 
of the director ^13^, a new set of elastic terms must be used [14^^15 : 

= ^L[{diajf + ^L^idittif + ^L'^aiaj{diak){djak) + ^i^4(e^J7ctt/c9^%•)^ 

where unlike the Frank elastic parameters the elastic constants L'- do not depend on the 
degree of order. To keep the number of elastic parameters minimal, we have retained among 
all possible terms quadratic in the derivative only those that are non-vanishing in the limit of 
a fixed degree of order. The total free energy, with included chiral coupling, is then derived 
in the form 

/ = -pC^^a^ + \pCa^ (2) 
+ ^p^Li{diajf + ^P^L2{diai)^ + ^P^L^ieijkakdiaj + ^o)^ (3) 



p1 



(4) 



x[p{p-Pc)r + i^L,{dip)\ (5) 
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where C is a positive Landau constant describing the isotropic-nematic phase transition, p* 
is the transition density, qo is the wave vector of the bulk cholesteric phase, and i^i,2,3 can 
be related to the Franck elastic constants, while x specify the rigidity of density 

variations. A more restrictive form of this free energy with no chiral interactions was derived 
in 
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The nonlinear density factor in the first term of the total free energy guarantees that the 
bulk nematic ordering stays limited to |a| < 1. In the ordering part of the free energy, ([2|, 
we have taken into account that each term is proportional to the number of molecules, i.e., to 
the local density p. The elastic free energy density ^ is however proportional to p^, as is the 
case for any interaction energy density. The continuity requirement ([T]) is taken into account 
by means of the penalty potential Q proportional to a coupling constant G. The optional 
density of chain beginnings and ends p^/^o will not be considered in this paper (see iq] 
also for a self-consistent distribution of p^). The density part ^ exhibits two singularities 
guaranteeing that the density stays positive and lower than the maximal packing density pc. 

Method of solution. - The equilibrium configuration of both constitutive fields, i.e., the 
density field p and the non-unit nematic director field a, are determined by minimizing the 
free energy at the constraint of a fixed global polymer mass, i.e. fixed polymer length. The 
minimization is performed by following a quasi-dynamic evolution of the director and density 
fields of the type described in 10 . 

We use a tangentially degenerate boundary condition for the director, i.e., the director 
is everywhere parallel to the surface of the sphere, while it is allowed to rotate freely in the 
tangential plane. For the density boundary condition at the confining shell we consider two 
separate cases: either no adhesion is assumed, with density set to zero at the inner surface 
of the shell, mimicking a short range repulsive interaction, or the density at the surface of 
the shell is fixed to be equal to the average density. This mimicks a short range attractive 
interaction when compared to the previous zero polymer density boundary condition. The 
initial condition is a homogeneous density field (except for the step at the boundary) and 
a(r) = plus a small random perturbation for the director field. The equations are solved 



by an open source finite volume solver 19 on a 50 x 50 x 50 cubic mesh (size of the box 
containing the sphere), which is adjusted near the surface of the sphere to define a smooth 
boundary. 

Units. - We use a reference nematic correlation length (the characteristic size of the 
defect core) as the length unit, defined as = yLopoJC^ where Lq and po are fixed reference 
elastic constant and polymer density, respectively. Length is thus expressed as r = and 
density as p = pop, where ~ denotes the dimensionless quantity. The correlation length of the 
nematic DNA was picked at 8 nm [Tg , which presents an upper bound for the bulk DNA 
experiments, and serves as the connection between the length scale of our simulation and 
the physical length scale. Expressing the free energy density ([2])-([5| in units of poC, the 
parameters of the model appear in the following dimensionless form, denoted by tilde: C = 1, 
Li = LiPo/C(,Q^ G = Gpo/C(,Q^ X = x/C'Po, Lp = Lppo/C^o- From now on all quantities 
will be dimensionless and ~ will be omitted. 

Results. - In what follows we present the steady state solutions of the quasi-dynamic 
evolution of the director and density fields corresponding to direct solutions of the Euler- 
Lagrange equations. We take the following dimensionless values for the parameters entering 
our model: Li = 1, L2 = 0, L4 = 1, x = 3 • 10~^, Lp = 0.1. The nematic transition threshold 
density is chosen p* = 0.5 and the tight packing density pc = 5. The coupling constant G was 
set to G = 3 and larger where required, so that the constraint ([T]) was satisfied and further 
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(a) qo =0 




(e) qo = 0.35 




(b) go = 0.1 




(f) go = 0.4 




(c) qo = 0.3 




(g) go = 0.6 




(h) qo = l 



Fig. 1 - (color online) Sequence of increasing chiral strengths qo of the nematic polymer in the sphere 
of dimensionless radius 32 (see main text) at an average density of p = 1. The density field is fixed 
to zero on the surface of the confining sphere. The I-N transition density is p* — 0.5, whereas pc — ^ 
is the tight packing density. The chirality go = 1 corresponds to the bulk cholesteric pitch of 50 nm 
chosen to accentuate the trends. The slight truncation of the torus in (a)-(b) is an artifact of the 
(rectangular) computational mesh. For clarity, only a representative part of the director field is shown 
in (e)-(g) and none is shown in (h). 



increase of G had no influence. Configurations in a sphere with radius of 32 are presented, 
corresponding roughly to physical radius of 250 nm for the chosen upper bound value of the 
nematic correlation length. 

We first study the density and the orientational order of a fully ordered case at increasing 
chiral interactions, Fig. [l] where the pitch of the bulk cholesteric phase, 27r/(7o, varies from 
oo (a) to 50 nm (h) . The polymer- wall interaction is assumed to be repulsive leading to zero 
surface density boundary condition (impenetrability). Moderate chirality obviously leads to 
a twisted toroidal orientational ordering, where the nematic director of the polymer circles 
around the centerline of the torus in the polar direction. The twist deformation increases 
with increasing the chiral strength, (a)-(b), with the morphology of the toroidal condensate 
remaining unaffected. As the chirality increases further, (c), the torus gets twisted, becoming 
folded [It] and more globule-like, while the director in the center is aligned with the symmetry 
axis of the original torus and is thus regular everywhere. Increasing the chiral coupling even 
further in (d) we observe the evolution of the twisted torus into a structure resembling a 
simple link. At extreme chiral strengths the condensate finally breaks up into a tube-like 
network filling space, (e)-(h). The director runs along the central line of the tubes and winds 
helically around it, while the tubes meet in a configuration that allows the director to be 
regular everywhere in the condensate. 

In order to see the details of the nematic transition of the polymer chain inside the spherical 
enclosure we now perform a density scan of the minimizing solutions, assuming that the 
enclosed polymer length and thus the average density p varies. Again we assume zero surface 
density boundary condition (impenetrability). Fig. [2] presents a density sequence for p = 



Daniel Svensek and Rudolf Podgornik: Polymer nematic condensation 



5 




(a) p = 0.35 (b) p = 0.4 (c) p=l (d) p = 3.5 



Fig. 2 - (color online) Sequence of increasing average densities p of the polymer in the same sphere 
as in Fig. [l] no chirality. The density field is fixed to zero on the surface of the confining sphere. The 
I-N transition density is p* = 0.5, whereas pc = 5 is the tight packing density. The contour denotes 
p = 2.5. The transition threshold is lowered and the polymer is fully condensed into the high density 
torus as soon as just above the threshold (b). 



0.35,0.4,1,3.5. One observes that the I-N transition leads to a breaking of the spherical 
symmetry of the density and the orient ational fields, yielding immediately a fully formed 
condensate in the form of a torus. This torus floats inside the sphere as long as it can, i.e., as 
long as the average density is small enough so that it does not touch the inner surface of the 
sphere. After that, it is pressed against the surface and deforms into a spheroid and finally 
into a sphere [5j[7|[8]. This sequence of events leading to different scalings of the external 
radius of the torus with its volume is similar to the case treated in ^18j within a different 
context. The short range steric interaction between the polymer condensate and the spherical 
enclosure preserves the symmetry of the aggregate at all packing densities. 

The situation is different with the boundary condition that sets the density at the enclosing 
surface equal to the average density, simulating an attractive short range surface interaction 
between the polymer condensate and the enclosing wall, leading to a non-vanishing surface 
density. This particular choice allows us to use the solving algorithms that are already working. 
First of all, we note in this case. Fig. |3j that the symmetry breaking at the FN transition is 
different then in the case of the polymer excluding surface. The condensate in fact wants to 
approach the wall and when its outer radius is smaller then the inner radius of the enclosure, 
it needs to break the polar symmetry of the condensed solution. Roughly, one could say 
that the condensate becomes glued onto the wall which deforms its original toroidal shape. 
At larger densities part of the symmetry of the toroidal condensate observed in the case of 
repulsive confining surface interactions is restored. On increase of the average density the 
interplay of volume elasticity of the polymer condensate, the FN interface, and the polymer 
interaction with the enclosing walls, leads to a complete restoration of the polar symmetry 
and the toroidal condensate becomes cup-like. Fig. |3](d). With this shape of the FN interface 
it can apparently minimize the total free energy subject to all the constraints. On increasing 
the average density even further, this cup-like toroid grows and eventually reaches the same 
spherical final state as in the case of a purely repulsive interaction with the confining wall, 
see Fig. [2](d). This of course makes sense since at large average density the condensate just 
tries to fill all the space available. The adhesion of symmetry broken shapes to the bounding 
surface has not been observed before in an approximate analysis of the encapsidated DNA 
toroids [5]. 

Moderate chirality does not have a pronounced effect on the transition but it does show 
up in the nematic director texture of the toroidal condensate. Fig. [4] shows the transition 
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(a) p = 0.5 (b) p = 0.6 (c) p = 0.75 (d) p = 2.5 



Fig. 3 - (color online) Sequence of increasing average densities p of the polymer in the same sphere as 
in Fig. [l] no chirality. The density is fixed to p on the surface of the sphere, mimicking an attractive 
polymer-surface interaction. The I-N transition density is p* — 0.5, whereas pc = 5 is the tight 
packing density. The contours denote p = 2.5, except in (c) where it denotes only a slightly larger 
value to make the plot sensible. For this boundary condition the I-N transition threshold is not 
lowered (a) as it is in Fig.|2] but the immediate condensation is there. The most pronounced effect 
of the boundary surface interactions is the breaking of the polar symmetry of the condensate. The 
toroidal condensate is adsorbed to the sphere and grows from there when the sphere is filled up. 



sequence for chirality qo = 0.1 and density fixed to p on the surface of the sphere. The 
contours show a winding helical configuration of the director around the circumference of the 
toroidal condensate. The exact position of the small incipient toroidal aggregate, Fig. [4](b), 
observed close to the I-N transition (where the torus is on its way to adsorb to the surface 
of the enclosing sphere), i.e. whether it is in close proximity to the boundary of the enclosing 
sphere or it floats somewhere within the sphere, is difficult to pinpoint exactly as the energy 
differences are very small and initial conditions of the calculation matter. Nevertheless, as 
the aggregate grows. Fig. [4](c)-(d), the final cup-like state becomes independent of the initial 
conditions and corresponds to a deep free energy minimum. Enhancing the chirality finally 
leads to a breakup of the condensate as shown in Fig. [l](f)-(h). 

Discussion. - The Landau-de Gennes theory of confined polymer nematic ordering pre- 
sented above describes the nematic ordering itself as well as the equilibrium shapes of the 




(a) p = 0.49 (b) p = 0.55 (c) p = 1.5 (d) p = 2.5 



Fig. 4 - (color online) Sequence of increasing average densities p of the polymer in the same sphere 
as in Fig. [l] with chirality = 0.1, corresponding to a bulk cholesteric pitch of 0.5 pun. The density 
is fixed to p on the surface of the sphere. The I-N transition density is p* = 0.5, whereas pc = 5 is 
the tight packing density. The contour denotes p = 2.5, except in (d) where it denotes only a slightly 
larger value to make the plot sensible. The effect of chirality is clearly seen in the texture of local 
director that now winds around the circumference. 
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ordered condensate for long chiral polymers with specific short range interactions with the 
enclosing wall. These interactions can be either repulsive, leading to polymer exclusion from 
the vicinity of the bounding surface and thus to vanishing polymer density at the spherical 
boundary, or effectively attractive with a corresponding finite boundary density. Within this 
approach we were first of all able to describe the effects of polymer chirality, which in the 
extreme case course-grain the nematic condensate into a tube-like network filling space, that 
then appears like a complicated arrangement of ordered domains as indeed seen in recent 
experiments 9 . In a weaker form chirality shows up in a winding helical configuration of the 
director of the toroidal condensate. Furthermore, in the case of effective attractive interac- 
tions between the condensate and the enclosing spherical shell, on formation the condensate 
needs to break the polar symmetry of the original toroidal shape in order to adsorb onto the 
enclosing wall. The broken polar symmetry of the original toroidal aggregate is later restored 
at higher average densities leading to a cup-like toroid that on increase of the average polymer 
density eventually reaches the same spherical final state as in the case of a purely repulsive 
interaction with the confining wall. 

The theory of nematic polymer condensation and ordering in a spherical confinement, de- 
scribing the DNA condensation within a virus capsid, presented above is a relevant alternative 
to computer simulations |20 , not demanding huge computational resources. It is applicable 
not solely to DNA condensation in bacteriophage capsids but to any nematic polymer-filled 
virus-like nano particles and should thus find a wide range of possible applications. 

* * * 
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